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9^. , ABSTRACT 

We study the slow neutron capture process (s process) in Asymptotic Giant 
Branch (AGB) stars using three different stellar evolutionary models computed 
for a 3 Mq and solar metallicity star. First we investigate the formation and 
the efficiency of the main neutron source: the ^^C(a,n)^^0 reaction that occurs 
in radiative conditions. A tiny region rich in ^'^C (the ^^C pocket) is created by 
proton captures on the abundant ^^C in the top layers of the He intershell, the 
zone between the H shell and the He shell. We parametrically vary the number of 
protons mixed from the envelope. For high local protons over ^^C number ratio, 
p/^^C > 0.3, most of the ^^C nuclei produced are further converted by proton 
capture to ^^N. Besides, ^"^N nuclei represent a major neutron poison. We find 
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that a linear relationship exists between the amount of ^^C in the He intershell 
and the maximum value of the time-integrated neutron flux. Then we generate 
detailed s-process calculations on the basis of stellar evolutionary models con- 
structed with three different codes, all of them self-consistently finding the third 
dredge up, although with different efficiency. One of the codes includes a mecha- 
nism at each convective boundary that simulates time-dependent hydrodynamic 
overshoot. This mechanism depends on a free parameter /, and results in partial 
mixing beyond convective boundaries, the most efficient third dredge up and the 
formation of the ^^C pocket. For the other two codes an identical ^^C pocket is 
introduced in the post-processing nucleosynthesis calculations. The models typi- 
cally produce enhancements of heavy elements of about two orders of magnitude 
in the He intershell and of up to one order of magnitude at the stellar surface, 
after dilution with the convective envelope, thus generally reproducing spectro- 
scopic observations. The results of the cases without overshoot arc remarkably 
similar, pointing out that the important uncertainty in s-process predictions is 
the ^^C pocket and not the intrinsic differences among different codes when no 
overshoot mechanism is included. The code including hydrodynamic overshoot 
at each convective boundary finds that the He intershell convective zone driven 
by the recurrent thermal instabilities of the He shell (thermal pulses) penetrates 
the CO core, producing a He intershell composition near to that observed in H- 
deficient central stars of planetary nebulae. As a result of this intershell dredge 
up the neutron fluxes have a higher efficiency, both during the intcrpulsc peri- 
ods and within thermal pulses. The s-element distribution is pushed toward the 
heavier s-process elements and large abundances of neutron-rich isotopes fed by 
branching points in the s-process path are produced. Several observational con- 
straints are better matched by the models without overshoot. Our study need 
to be extended to different masses and metallicities and in the space of the free 
overshoot parameter /. 

Subject headings: nucleosynthesis, abundances — stars: AGB, evolution 



1. Introduction 

When helium is exhausted in the center of low- to intermediate-mass stars (M ^ 8 
Mq) the core becomes degenerate (Paczyhski 1970) and burning processes do not occur 
there anymore. Instead, H and He burn alternately in shells while the star ascends the 
Asymptotic Giant Branch (AGB). In this phase the stellar structure consists of, from the 
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center outward, a degenerate CO core, the He-burning shell, a thin (10~^ — 10"^ M©) zone 
between the H shell and the He shell (hereafter He intershell), the H-burning shell, and a 
large convective envelope, which suffers from strong stellar winds. When the envelope mass 
is reduced to less than ~ fO^^ Mq the star leaves the AGB, a Planetary Nebula is formed 
and the CO core cools towards the White Dwarf phase. In 10-20% of cases the residual H- 
rich envelope is completely removed by processes associated with a late He-shell instability 
during the pre- White Dwarf phase (Herwig et al. 1999), leading to the formation of the 
H-deficient central stars of planetary nebulae. 

During the AGB the H shell dominates the energy production for most of the time. 
Hydrogen is transformed into He, the He intershell consequently grows and is progressively 
compressed and heated until the temperature and density become high enough that He 
burning is triggered in the bottom layers. The thermonuclear runaway, also known as thermal 
instability, or thermal pulse, (Schwarzschild & Harm 1965) generated by this sudden release 
of energy causes almost all the He intershell to become convective (we refer to this as the pulse 
driven convective zone, PDCZ) while the envelope expands and the H shell cools. Within 
the PDCZ partial He burning produces large amounts of carbon. The PDCZ quenches 
after a time of a few ten to a few hundred years, while He burning continues radiatively 
for another few thousand years, during which the H shell is inactive. Then, the envelope 
contracts and shell H burning starts again. The cycle is repeated for a few up to possibly 
100 times, the total number of thermal pulse depending on the initial stellar mass and 
the mass loss rate, with thermal periods of a few lO^uptolO^ yr. Recent models of the 
AGB phase have been computed among others by HoUowell Sz Iben (1988); Boothroyd & 
Sackmann (1988a,b,c,d); Lattanzio (1986, 1987, 1989); Vassihadis & Wood (1993); Blocker 
(1995); Forestini & Charbonnel (1997); Straniero et al. (1997); Wagenhuber & Groenewegen 
(1998); Mazzitelh, D'Antona, & Ventura (1999); Mowlavi (1999); Herwig (2000). 

As illustrated in Fig. 1, after a limited number of thermal pulses, when the mass of the 
H-exhausted core is above ~ 0.6 Mq, the convective envelope penetrates in the top layers of 
the He intershell bringing to the surface newly synthesized He, C and elements produced by 
neutron captures. This recurrent phenomenon is called third dredge up (TDU). The TDU 
is responsible for the carbon enrichment shown by MS, S and C(N) stars. Unfortunately, 
the reproducibihty of the TDU phenomenon is a main problem with AGB computations, 
connected to the more general astrophysical problem of the treatment of turbulent convection 
in stellar interiors. Different evolutionary codes do not reproduce the same results: some 
codes found a large amount of TDU at solar metaUicity only for intermediate-mass stars 
(5 Mq < M < 8 Mq, Iben 1975a), others found TDU to occur in low-mass stars (M < 
5 Mq), but only for low metallicities: 1/200 of solar (Iben & Renzini 1982a) and 1/20 of 
solar (Boothroyd & Sackmann 1988d). Lattanzio (1989) and Straniero et al. (1997) found 
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TDU to occur also in stars with initial mass as low as 1.5 Mq and solar metallicity. The 
amount of TDU is also connected to the evolution of the envelope depending, for example, 
on the temperature at its base (Wood 1981) and its total mass (Straniero et al. 1997). The 
inclusion of extra mixing in the computation of TDU, however, can affect these results. The 
lower limit of 1.5 Mq at solar metallicity for TDU to occur is consistent with inferences 
derived from observational constraints of carbon stars (Groenewegen, van den Hoek, & de 
Jong 1995). Frost & Lattanzio (1996) studied the effect of the numerical treatment on the 
occurrence and magnitude of TDU, while Mowlavi (1999) tested the effect of using an extra 
mixing mechanism at the base of the convective envelope and found that TDU occurs in 
his models only if this procedure is apphed. Recently, Pols & Tout (2001) showed that 
when the AGB evolution is computed by solving simultaneously the structure equations and 
the diffusion equations for changes in chemical composition the TDU strongly depends on 
numerical choices, such as how the diffusion coefficients are averaged. 

Another debated problem is a satisfactory description of mass loss during the AGB 
phase. It is believed that the main mechanism responsible for mass loss in red giants is 
radiation pressure on dust grains. Detailed calculations in conjunction with new observations 
will address this problem. Most recent works include Wachter et al. (2002) and Olofsson et 
al. (2002). In the past an often employed way to include mass loss has been to use Reimers 
(1975) 's formula, in which the mass loss rate depends on the luminosity, radius and mass 
of the star and is proportional to a free parameter, rj, whose value can vary (typically from 
~ 0.3 to ten). Other semi-empirical relations for the mass-loss rate have been proposed, in 
particular by Vassihadis & Wood (1993); Blocker (1995) and Arndt, Fleischer, & Sedlmayr 
(1997). Marigo (2002) investigated the effects of molecular opacities on the evolution of AGB 
stars as the C/0 ratio grows from below to above unity and found that in the carbon-rich 
models a notable cooling results and hence an early onset of the superwind is expected in 
the framework of the Vassihadis & Wood (1993) 's formulation of the mass loss rate. 



1.1. The s process 

The He intershell is the site of the s process, the neutron capture nucleosynthesis that 
occurs when neutron number densities Nn are of the order of 10® — 10^^ n cm~^. In this range 
of Nn when an unstable nucleus is produced along the s-process chain of neutron captures, 
it typically decays rather than capture another neutron (Burbidgc et al. 1957). However, 
for relatively high values of neutron density or temperature several branchings along the 
s-process path are open and neutron- rich nuclei can be produced (Ward et al. 1976). As 
a "rule of thumb" the s process distribution follows the law a{A)Ns{A) — constant, where 
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(t{A) is the neutron capture cross section of nucleus A and Ns{A) its s-process abundance 
by number. This rule is only valid locally because of the presence of the neutron magic 
nuclei at iV = 50, 82, 126. These nuclei have very low neutron capture cross sections, 
which make them behave as bottlenecks for the neutron fiux. As a consequence steps in the 
otherwise a{A)Ns{A) = constant distribution are found in correspondence to neutron magic 
nuclei. The depth of these steps is a function of the total time-integrated neutron exposure 
r = Jq NnVxdt (where vt is the thermal velocity), so that different neutron exposures lead 
to different a{A)Ns{A) distributions (see Clayton 1968, page 568). In AGB stars heavy s- 
process elements are produced by the s process and then mixed to the stellar surface by TDU 
where they have been observed (for general references see Smith & Lambert 1990; Wallerstein 
et al. 1997; Wallerstein & Knapp 1998; Busso et al. 2001; Abia ct al. 2001). In particular, 
AGB stars are responsible for the s-process main component (Kappeler et al. 1989; Arlandini 
et al. 1999; Travaglio ct al. 1999), i.e. a significant fraction of the Galactic production of 
elements from Sr to Pb, as well as for the strong component (Clayton & Rassbach 1967) 
feeding a large fraction of solar ^osp]^ (Gallino et al. 1998; Goriely & Mowlavi 2000; Travaglio 
et al. 2001). 

The ^^C(«,n)^^0 reaction is activated at low temperatures (^ 0.8 x 10^ K), and plays 
the major role as neutron source in AGB stars (HoUowell & Iben 1988; Gallino et al. 1988; 
Kappeler et al. 1990). However, a higher amount of ^^C than present in the H- burning 
ashes is needed to reproduce the observed enhancements of heavy elements. To enable the 
formation of a region rich in ^^C, the ^^C pocket, a few protons from the envelope must 
penetrate into the He intershell and then react with the ^^C nuclei abundantly present there. 
The possibility of penetration of protons from the envelope directly into the convcctive PDCZ 
was found by Ibcn (1976) to be severely inhibited by an entropy barrier induced by the still 
active H-burning shell. A favorable location for the mixing of protons is the sharp H/He 
discontinuity that is left over after the occurrence of TDU while the H shell is temporarily 
extinguished (Iben & Renzini 1982b). Recently, models including time- dependent overshoot 
(Herwig et al. 1997) motivated by hydrodynamic simulations by Freytag, Ludwig, & Steffen 
(1996) (see also Singh, Roxburgh, & Chan 1998), and models with rotation (Langer et 
al. 1999) have been able to create a proton rich layer with Mpocket ~ 10~^ — 10~^ Mq at 
the top of the He intershell after the end of the TDU. A similar effect was obtained by 
imposing a velocity gradient profile beyond the border of the convective envelope (Cristallo 
et al. 2001). Another description of turbulent perturbation of stable layers below convective 
zones can be performed in the framework of internal gravity waves (Denissenkov & Tout 
2002). In summary, the occurrence of hydrogen penetration at the top of the He intershell 
seems plausible, but the mass involved and the resulting ^^C abundance profile are still to 
be treated as relatively free parameters, constrained by the rates of proton captures on the 
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carbon isotopes. 

The s- process mechanism in AGB stars is visuahzed in Fig. 1. After less than a few- 
thousand years ^^C and ^^N are produced by proton capture on the abundant ^^C in the top 
layers of the He intershell. Subsequently, after about 10,000 yr, the ^^C(Q;,n)^^0 reaction is 
activated during the interpulse phase in radiative conditions at a relatively low temperature: 
~ 0.9 X 10^ K (Stranicro et al. 1995). Before the end of the interpulse period all ^^C has burnt 
and the pocket has become enriched in s-processed material. This s-process nucleosynthesis 
occurs in radiative conditions, hence there is no interaction among the different layers of 
the pocket. Each zone can be treated separately and the efficiency of the neutron flux is a 
local property that only depends on the initial p/^^C ratio in that layer. The neutron flux 
lasts typically 20,000 yr and produces locally high neutron exposures, of the order of 10~^ 
mbarn"^. On the other hand, because the timescale is long, the neutron density remains 
low, reaching about 10^ n cm~^. At the end of the interpulse the pocket is engulfed by 
the following PDCZ and mixed with ashes from the H-burning shell and material from the 
previous PDCZ. A large amount of ^^Ne is present in the PDCZ as a product of the chain 
^^N(a,7)^^F(/9"'"z/)^^0 (a,7)^^Ne starting on the abundant ^^N from the H-burning ashes. In 
the convective zone the ^^Ne(a,n)^^Mg reaction can be marginally activated and a second 
neutron flux occurs. This second neutron burst is opposite in features to that in the pocket: 
it occurs on a timescale of a few years and it produces low neutron exposures (of the order of 
10~^ mbarn"^) with a high-peaked neutron density, up to 10^° n cm~^. This neutron burst 
does not typically contribute much to the overall production of s elements, however it affects 
the final abundances of isotopes connected to branching points. For an extensive review on 
the s process in AGB stars see Busso, Gallino, & Wasserburg (1999). 

We want to analyze the dependence of s-process model predictions on the stellar evolu- 
tionary code used to compute them. We compare the results of post-processing models based 
on three evolutionary sequences computed with three different codes, one of them including 
an overshoot mechanism. The sequences are different in the details of the description of phys- 
ical processes like mass loss, and of the numerical solution scheme adopted. To our purposes 
these computations represent a typical set of models of the same star constructed by different 
groups. The paper is organized as follows: in §2 we describe the stellar evolutionary codes 
and the evolutionary sequences we used to compute the nucleosynthesis. In §3 we describe 
the nucleosynthesis codes, the formation and efficiency of the ^^C neutron source, and the 
neutron fluxes and final s-process distributions as computed with the different evolutionary 
sequences. In §4 we compare the final results with a number of observational constraints. 
In §5 we outline our conclusions. A preliminary study was presented by Lugaro & Herwig 
(2001). 
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2. The evolutionciry codes 

To establish the main effects and dependencies we consider only the case of a star of 
initial mass 3 Mq and solar metallicity. The most recent s-process studies (Gallino et al. 
1998; Busso et al. 1999; Goriely & Mowlavi 2000) have shown that there is an important 
variation in the resulting s-process distribution with stellar metallicity. The ^^C in the 
pocket is a primary neutron source having formed from H and the ^^C resulting from He 
burning, hence for decreasing metallicities heavier and heavier elements are produced because 
there are more neutrons per iron seed. This property has been shown to have important 
consequences in the study of the s process at different metallicities (Busso et al. 2001) and 
hence on the galactic chemical evolution of heavy elements (Travaglio et al. 1999). We will 
not discuss this point to a large extent since we only consider models of solar metallicity. 
Busso et al. (2001) demonstrated that a spread in the efficiency of the neutron release in 
the ^^C pocket is necessary to cover the spectroscopic observations at any given metallicity. 
Again, we will not discuss this result in detail because for sake of a simple comparison we 
perform s-process calculations using a given ^^C pocket, chosen as a relatively free parameter 
in the evolutionary sequences without overshoot, or as resulted from the computation in the 
evolutionary sequence with overshoot included. 

We consider three different stellar evolutionary codes: the Prascati RAphson Newton 
Evolutionary Code (FRANEC, Straniero et al. 1997), the Mount Stromlo Stellar Struc- 
ture Program (MSSSP, Wood & Zarro 1981; Frost & Lattanzio 1996), and the EVOL code 
(Blocker 1995; Herwig 2000). Since the codes are independent they contain several differ- 
ences, e.g. with respect to the numerical treatment, opacities and nuclear reaction rates. The 
most important difference to be noted here, however, is the treatment of convective insta- 
bilities and the associated mixing, which influences the calculation of the TDU phenomenon 
and the formation of the ^^C pocket. 

In the framework of the mixing length theory mixing velocities should drop to zero at 
the boundary of a convective layer where the radiative temperature gradient Vrad equals the 
adiabatic temperature gradient V ad-, ^rad ~ ad = 0, defining the Schwarzschild boundary. 
In reality, because of inertia, convective eddies are more likely to have non-zero velocity and 
to be dragged across the Schwarzschild boundary. This idea is the basis of the concept of 
mechanical overshoot. However, when convection penetrates into a chemically inhomoge- 
neous region a discontinuity in the opacity forms at the stability border. When the opacity 
in the convective region is larger than the opacity in the radiative region then Vrod ~ Vad > 
at the convective border seen from the radiative zone. In this condition extra mixing is in- 
duced by the change of composition and this causes the border to advance into the radiative 
zone. A famous example of such situation is found at the outer border of the convective core 
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during central He-burning: the opacity increases, due to the conversion of He into C and O, 
and the convective region grows in mass (Paczyhski 1970; Castellani, Giannone, & Rcnzini 
1971). The same phenomenon is encountered during the TDU phases (Becker & Iben 1979; 
Castellani, Chicfii, & Straniero 1990): the H-rich convective envelope penetrates into the 
He-rich intershell layer that has a substantially lower opacity. 

To handle this special condition in AGB models the MSSSP code uses a numerical 
scheme which may involve mixing of an additional stable mass shell after each structure 
iteration: "Suppose mesh point j is the outer (inner) edge of a formally convective zone. 
By evaluating the radiative to adiabatic ratio Vrad/^ad at mesh point j — 1 (j -|- 1), we can 
extrapolate to find {'Vrad/'Vad)extrap at j + 1 (j — 1). If this value is greater than unity, then 
we incorporate the j + 1 (j — 1) mass shell into the convective zone, even though the shell may 
have Vrad < "^ad when calculated formally. This procedure is executed, at each convective 
boundary, each time the abundances are mixed over a convective region." (Lattanzio 1986). 
Note that in the MSSSP code convective regions are mixed after each structure iteration. A 
further modification to MSSSP was described by Frost & Lattanzio (1996) to allow for the 
case where the increasing depth of the convective envelope (with each iteration) was followed 
by a retreat of the envelope convection, leading to failure to converge. When this occurs the 
iterations are treated as capturing the essential time-dependent response of the structure to 
the induced mixing. Hence the greatest depth of the convective envelope is taken as the 
mixed region, which cannot later be unmixed during subsequent iterations for this model. 
This mixed abundance is taken as the composition for subsequent iterations, perhaps leading 
to convergence with a shallower convective zone but with the deeper mixed region. 

No special algorithms are adopted by the FRANEC code version used here to handle 
the discontinuity in the chemical composition and the consequent jump of the opacity that 
forms at the base of the convective envelope. Since both MSSSP and FRANEC assume 
instantaneous (homogeneous) mixing, neither of them can find the proton mixing into the 
He intershell at the end of TDU needed for the formation of a ^^C pocket. The main difference 
between the two codes is the amount of TDU mass, which is higher in the MSSSP case. 

In EVOL a time- and depth-dependent overshoot mechanism is included at the convec- 
tive borders within the star, which induces partial mixing beyond the boundary of homoge- 
neously mixed regions. This scheme mimics mixing due to an exponential decay in space of 
the convective velocity field into the stable layers, as observed in hydrodynamic simulations, 
following the prescription of Freytag et al. (1996) who performed two-dimensional hydro- 
dynamic simulations of time-dependent compressible convection in order to study shallow 
stellar surface convection zones. A similar behavior of the decay of the convective velocity 
field outside the convection zone has been found by Asida & Arnett (2000) in two-dimensional 
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hydrodynamic simulations of the oxygen-burning shell in massive stars. The scale height of 
these type of mixing remains a free parameter (/) for the convective boundaries considered 
here. This type of hydrodynamical overshoot has also been used by Mazzitelli et al. (1999) 
and Schlattl & Weiss (1999). 

Herwig (2000) showed that this algorithm, when applied to thermally pulsing AGB 
models, is relevant to address some major observational properties. Time-dependent hydro- 
dynamic overshoot is a process that naturally leads to the formation of a ^^C pocket needed 
for the occurrence of the s process over a range of 10^^ — 10^^ Mq. It increases the TDU 
efficiency and therefore may better account for low-luminosity carbon stars (Herwig, Blocker, 
& Driebe 2000). It leads to intersliell dredge up of C-0 core material in the helium intershell. 
This enhances the C and O intershell abundances as required by models such as those pre- 
sented by Herwig et al. (1999) that reproduce observations of hot H-deficient post- AGB stars 
of type PG 1159 (Dreizler & Heber 1998) and central stars of planetary nebulae of spectral 
type [WC] (Leuenhagen & Hamann 1998). All these effects depend on the value of a free 
overshoot parameter /. Here we want to check whether the inclusion of hydrodynamic over- 
shoot is in contradiction with other observational informations, in particular those connected 
with the activation of the neutron sources and the consequent s-process nucleosynthesis. We 
use a value of /=0.016 at every convective boundary, except at the base of the convective 
envelope where we use /=0.128 for the computation of the ^^C pocket. The original value of 
0.016 was chosen to model main sequence core overshoot in order to reproduce the observed 
width of the main sequence (Schaller et al. 1992). However, with this choice the mass of the 
^^C pocket and consequently the total mass of s-process material is too small to account for 
observations (as already pointed out by Goriely & Mowlavi 2000). 

Cristallo et al. (2001) also find partial mixing of protons below the base of the convective 
envelope with the FRANEC code. They achieve this by means of a time-dependent mixing 
scheme (see Chieffi et al. 2001) with mixing velocities derived in the framework of the mixing 
length theory. To handle the Vrad — Vad discontinuity condition during TDU Cristallo et 
al. (2001) assume that the mixing velocity drops to zero below the Schwarzschild boundary 
of the convective envelope with an exponential decline in space. This choice induces both a 
deeper TDU and the formation of a tail of the internal H profile, thus allowing the subsequent 
formation of the ^^C pocket. Also in this case, the efficiency of the TDU and the extension 
of the resulting ^^C pocket depend on the efficiency of the assumed exponential dechne of 
the velocity field, i.e. a free overshoot parameter fy. The main difference with respect to 
the scheme adopted by EVOL is that in Cristallo et al. (2001) the overshoot is induced 
only where there is a Vrad — ^ad discontinuity, i.e. at the base of the convective envelope 
during TDU. We will also briefly comment on the s-process results obtained by using as 
inputs the ^^C pocket and the stellar features of the FRANEC stellar model of 3 Mq and 
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solar metallicity with Reimers (1975) 's parameter rj — 1.5 and overshoot parameter /„ = 
0.1 presented by Cristallo et al. (2001). This computation was followed for 17 pulses with 
TDUs, resulting in a final envelope mass of 1.3 Mq. Because of the algorithm of overshoot 
included at the base of the convective envelope, the dredge-up mass is higher than in the 
FRANEC standard case, reaching up a total of 0.06 Mq and C/0 is above unity in the last 
5 pulses. The ^^C pocket was found to form over a mass of about 4x10"^ Mq. 

Mass loss has also been modeled in different ways: the FRANEC code uses Reimers 
(1975)'s prescription with rj = 1.5 throughout the whole evolution, the MSSSP evolution 
has been computed with Reimers (1975) 's mass loss with 1] = 1 until the beginning of the 
AGB phase and afterwards using the formula proposed by Vassiliadis & Wood (1993), which 
accounts for a superwind phase. The EVOL code uses Reimers (1975) 's mass loss with rj — 1 
from the pre-Main Sequence to the point where C/0—1 in the envelope. Afterwards it uses 
the Arndt et al. (1997) rate for carbon stars. 



2.1. The stelleir evolutionciry sequences 

We computed s-process detailed nucleosynthesis for a stellar model of 3 Mq and solar 
metallicity making use of the structural, thermodynamical and chemical characteristics ob- 
tained from the three different evolutionary codes described above. The FRANEC sequence 
ends because of convergence problems after 25 thermal pulses with TDU, when TDU is van- 
ishing because of the low envelope mass, the MSSSP sequence also stops for convergence 
problems after 18 pulses with TDUs, with an envelope mass similar to that of the FRANEC 
case. For the EVOL sequence the adopted mass loss law results in a low mass loss rate. Af- 
ter reaching the 13*'^ TDU we stopped the sequence under the assumption that an efficient 
superwind phase will remove the remaining envelope within a few thermal pulses. 

In Fig. 2 we present some quantities related to the evolution of the convective envelope 
during the AGB phase as functions of the mass of the H-exhausted core Mh: the mass 
dredged up by each TDU episode, the cumulative mass dredged up by all previous TDU 
episodes, the mass of the envelope and the C/0 ratio in the envelope. Each point corresponds 
to a thermal pulse followed by TDU. The dredged- up mass and the envelope mass are 
determined by the different choices made in each code described above. Because of the 
different mixing algorithm and assumption on overshoot, the amount of dredged up mass 
grows from FRANEC, which does not include any treatment of the discontinuity of the 
adiabatic to radiative gradient ratio at the base of the convective envelope, to MSSSP, 
which attempts to find the neutrality in that ratio, to EVOL, which includes an algorithm 
for hydrodynamic overshoot. For the three models, TDU starts after the 9*^'^ lO*'^' and P* 
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thermal pulse respectively, when the core mass has reached ~ 0.61 — 0.63 Mq. The AGB 
sequence with EVOL starts with a larger core mass because of overshoot during the main 
sequence. In the 3 Mq model presented by Hcrwig (2000) the third dredge up starts at 
the 4*^^ pulse. The difference is due to the different value of / applied at the base of the 
convective envelope. An interesting consequence of the deep TDU experienced by the EVOL 
model is that after the first six pulses the core mass does not grow anymore, rather decreases 
with time (Herwig, Schonberner, & Blocker 1998). 

The evolution of the envelope mass is determined by the choice of the mass-loss law. For 
FRANEC, the Reimers (1975)'s formula determines a continuous decrease of the envelope 
mass, for MSSSP, which uses the Vassiliadis & Wood (1993)'s mass- loss law, the mass of 
the envelope is kept quite constant until a superwind sets in during the last two pulses. In 
EVOL, which uses Reimers (1975) 's for the first 6 pulses and then Arndt et al. (1997) when 
the envelope is carbon rich, the mass loss is not as large as in the FRANEC and MSSSP 
cases and the envelope mass is about double than in the other two cases right to the end of 
the computed evolution. The C/0 ratio at each pulse is determined by the mass dredged up 
and the envelope mass, and also by the C and O abundances in the intershell which are both 
larger in the EVOL case because of intershell dredge up. All the models become carbon rich 
in the envelope after a given number of dredge up events. 

In Fig. 3 we compare some interesting features of the He intershell as functions of the 
core mass: the interpulse period, the extent in mass of the intershell convective region, the 
overlap factor r between subsequent pulses and the bottom temperature of the pulse driven 
convection zone. On this last quantity depends the activation of the ^^Ne neutron source. 
While in the FRANEC and MSSSP cases the interpulse period decreases with time, and is 
very similar to previous AGB models (such as those presented by Boothroyd & Sackmann 
1988c), in the EVOL case the interpulse period increases with time: from ~ 60,000 yr to ~ 
80,000 yr. This difference is a consequence of fact that in EVOL the core mass decreases 
with time, since the lower the core mass, the larger the interpulse period (Paczynski 1974). 
For the three stellar evolutionary codes examined, the mass of the pulse driven convection 
zone decreases with time and is higher in the EVOL case because the convection extends 
down in the CO core. The overlap factor r decreases with time in all cases. The decrease is 
steeper for cases with larger TDU. The maximum temperature at the base of the pulse driven 
convection zone grows with the pulse number in all three cases. The FRANEC trends are 
very close to those of Boothroyd & Sackmann (1988c). In EVOL the maximum temperature 
reached is 3.20 xlO^ K, with respect to the maximum value of 3.05 xlO^ K reached in the 
FRANEC and MSSSP case. This higher temperature is a consequence of the hydrodynamical 
overshoot applied at the base of the intershell convection zone. 
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3. The nucleosynthesis calculations 

We have used two different post-processing nucleosynthesis codes to different aims. To 
investigate the formation of the ^'^C podcet and evaluate the efficiency of the ^^C main 
neutron source we have made use of the Monash nucleosynthesis code; to perform detailed 
s-process calculations we have used the s-process code by Gallino et al. (f998). 

The Monash nucleosynthesis code reads a file of stellar structure inputs such as temper- 
ature and density that is produced by the MSSSP evolutionary code. The nuclear network 
is specifically designed for the detailed study of light-element nucleosynthesis in AGB stars. 
In the nucleosynthesis network there are 74 nuclear species: from neutrons and protons up 
to sulphur there are 59 nuclei, with another 14 iron group species to allow neutron capture 
on iron seeds. There is also an additional "particle" g which has the function of counting 
the number of neutron captures occurring beyond ^'^Ni hence giving an estimation of the 
resulting heavier elements s-process distribution. All proton, a, neutron captures and f3 
decays involving the species listed above are included in the nuclear network summing up 
to 506 reactions. The bulk of reaction rates are from the Reaclib Data Tables of nuclear 
reaction rates based on the 1991 updated version of the compilation by Thielemann et al. 
(1987). They have been updated to the latest new experimental or theoretical estimates 
(Lugaro 1998). The reaction network is terminated by a neutron capture on ^^Ni followed 
by an ad hoc decay with A = 1 s~^ producing the particle represented by the symbol g: 
^''■Ni(n,7)^^Ni ^^Ni + g. The Q-value assigned to this reaction is the value of the reverse 
reaction ^^Ni(7,n)^^Ni. Following the method of Jorissen & Arnould (1989) neutron captures 
on the missing nuclides are modeled by neutron sinks, meaning that the ^''S(n,7)''^^S and the 
^-^Ni(n,7)^^Ni reactions are given some averaged cross section values in order to represent 
all nuclei from ^^S to ^^Mn and from ^^Ni to ^°^Bi respectively. We obtained (J {^"^S) sink = 
21.1 mbarn at 8 keV (the typical temperature of ^^C burning) using a solar distribution 
for the abundances of nuclei from ^^S to ^^Mn. To compute a{^^'!^i)sink we used different 
abundance distributions of heavy elements computed in different phases of the s process in 
AGB stars of 3 Mq and solar metaUicity from the models of Gallino et al. (1998). We 
found that the results do not vary much: changing by up to 15% for different distributions. 
Therefore, we decided to use as typical the average of such values: (7{^^Ni)sink = 135 mbarn 
at 8 kcV. For stellar models of different masses or metallicity this parameter might have to 
be adjusted. The temporal variation of the neutron density is determined in each layer by 
the local equilibrium between neutron production and destruction. 

The s-process code follows a complete set of neutron captures on 404 nuclei from "^He 
up to Po, where the s-process path ends because this element is unstable against a decay. 
a-Capture reactions are considered on light nuclei up to Mg. This code uses structure and 
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thermodynamic inputs from computed stellar sequences. Reaction rates for neutron produc- 
tion have been taken, as in the Monash nucleosynthesis code, from Denker et al. (1995) for 
the ^^C(Q;,n)^^0 reaction and from Kappeler et al. (1994) for the ^^Ne(Q;,n)^^Mg reaction, 
excluding the contribution of the elusive resonance at 633 keV and using the lower experi- 
mental limit for the resonance at 828 keV. Rates for (n,7) reactions are updated to the latest 
estimates (Bao et al. 2000; Reifarth et al. 2001; Mutti et al. 2001). Rates for {n,a) and (n,p) 
reactions are taken from the same sources as in the Monash nucleosynthesis code. Decay 
rates are mainly from Takahashi & Yokoi (1987). 



3.1. The formation and the evolution of the ^^C pocket 

With the Monash nucleosynthesis code we have studied the formation of the ^^C pocket 
and the activation of the ^^C(Q;,n)^^0 neutron source reaction during a typical interpulse 
period. This is done to understand the effect of the ^^N(n,p)^''C neutron poison reaction 
and the subsequent recycling of protons by the ^^C(p,7)^^C reaction, possibly inhibited by 
the ^^C(p,7)^^N reaction and proton captures on other isotopes such as ^^O. We perform 
the calculations with the Monash nucleosynthesis code because it includes proton captures 
while the s-process code does not. We derive the main general features of the neutron flux 
in the ^^C pocket and their dependence on various parameters. 

We generated the ^^C pocket by introducing in the post-processing Monash nucleosyn- 
thesis code computation a certain number of protons below the H/He discontinuity at the end 
of the 10**^ TDU episode, which follows the 19*'^ thermal pulse in the 3 Mq model computed 
with MSSSR The initial mass fraction of ^^C in the He intershell is Xci2^=0.23, which is a 
very typical value for an advanced PDCZ (see e.g. Fig. 9 of Boothroyd & Sackmann 1988c). 
No ^^C or ^^N from H-burning ashes are present initially in the He intershell since they have 
been destroyed by a captures during the previous PDCZ. We performed one-zone model 
calculations by introducing a constant mass fraction of protons X*™* for each simulation in 
a region of about 10"^ Mq. The typical temperature of the top layer of the He intershell 
grows from about 2x10^ K to about 4.5x10^ K in the ~ 2000 yr following the end of the 
TDU. At 4.5 xlO^ K the lifetime of protons with respect to the ^^C(p,7)^^N reaction is < 
200 yr, so that all protons are captured before the onset of H-shcll burning, which occurs 
~ 3000 yr after the end of TDU. Carbon-13 is formed by proton capture on ^^C and ^"^N is 
produced by proton captures on the newly formed ^^C. The number of protons available 



^We indicate the mass fraction X of isotope ^^C with the notation Xci2, and equivalently for aU the 
other isotopes. 
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determines the ^^C and ■'^^N final concentrations. For the ■^^C(p,7)^^N reaction we used the 
rate from Caughlan & Fowler (1988), for the rate of the ^^C(p,7)^^N reaction we used the 
new evaluation from King et al. (1994), which in this range of temperature is 1.29 times 
higher than the previous evaluation by Caughlan & Fowler (1988). 

The results of these computations are presented in Fig. 4. As first shown by HoUowell 

& Iben (1988) and discussed in detail by Jorissen & Arnould (1989) and Goricly & Mowlavi 
(2000), for p/^^C < 1/10, the number of ^^C nuclei produced is essentially determined by the 
number of initial protons. Then, for increasing number of protons an increasing fraction of 
the ^^C nuclei suffer a second proton capture and are transmuted to ^"^N, until, at p/^^C > 1, 
a ^^N pocket emerges next to the ^'^C pocket. For X^^^ > 0.02, p/^^C > 1 and ^^C is strongly 
depleted, being transformed into ^^N. Since ^^N is a strong neutron absorber, the region of 
the pocket which is more efficient for the production of s-process nuclei is where ^^C is more 
abundant than ^^N. The results are affected by the choice of the ^^C(p,7)^^N rate; when the 
lower rate from Caughlan & Fowler (1988) is used the final amount of ^'^C is higher. For 
example, at Xp"** = 0.003, with the older rate we obtain 12% more ^^C and 21% less ^'^N . 

The ^^C(Q;,n)^^0 reaction occurs during the interpulse period because the temperature 
in the ^^C-pocket layers reaches ~ 10^ K before the start of the next pulse. The lifetime 
of ^^C with respect to the ^^C(a;,n)^^0 reaction is about 400 yr at this temperature. Thus 
^^C is totally destroyed before the onset of the next thermal pulse. With the MSSSP code 
the temperature in the ^^C pocket reaches 10^ K after about 43,000 yr (with p ~ 13, 000 g 
cm~^ and Xhc4 = 0.72) from the beginning of the interpulse period. The whole interpulse 
period lasts for about 52,000 yr. Thus, we concur with Straniero et al. (1995, 1997) that all 
the ^^C present in the pocket burns in radiative conditions. (However, it has to be noted that 
this might not occur during the earliest interpulse periods, see e.g. Fig. 4 of Herwig 2000). 
The resulting neutron exposure is shown in Fig. 4. A most important point is the role played 
by the ^^N(n,p)^''C reaction and the subsequent partial proton recycling (see also Gallino et 
al. 1998, and references therein). Nitrogen-14 nuclei indeed act as a strong neutron poison 
during the s process because of the relatively high cross sections of the ^^N(n,p)^^C reaction: 
recent measurements give a Maxwellian average cross section of 2.04±0.16 mbarn at 24.5 
keV (Gledenov et al. 1995) with a very nearly 1/v shape. Neutron capture cross sections 
of all the other CNO nuclei are typically smaller by about two orders of magnitude. The 
protons produced by this reaction are recaptured either by the abundant ^^C, and ^^C is 
produced again, or by the ^^C itself, which has a proton capture rate about ten times larger 
than that of ^^C. The overall effect can be considered as a partial recychng of neutrons 
together with a leaking of the neutron source. In no case is the proton recycling mechanism 
efficient enough to leave any ^^G present in the pocket at the end of interpulse period. When 
the mass fraction of ^"^N is smaller than that of ^^C, few neutrons are captured by ^^N and 
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the s-process efficiency is high. In this situation the small recychng effect produces slightly 
larger final neutron exposures (e.g. 10% higher for initial proton number of 0.001) than in 
the test case when proton captures on ^^C and ^"^C are switched off. When the mass fraction 
of ^^N is higher than that of ^'^C, most neutrons arc captured by ^^N to be transmuted into 
protons. In this situation the recycling effect produces larger final neutron exposures than if 
proton captures on ^^C and ^^C were not included (e.g. 30% more for initial proton number 
of 0.02). However, the s-process efficiency is low. For X*™* > 0.02 the abundance of ^^C, 
and hence the recycling effect, is much reduced and the neutron flux is unimportant. This 
explains why in Fig. 4 the total neutron exposure in the ^^C pocket reaches a maximum 
and then decreases as a function of Xp™*. In particular, we found that a maximum r of 0.43 
mbarn"^ can be achieved, for solar metallicity stars, when Xp™* is 0.003 (as also found by 
Goricly & Mowlavi 2000). In view of this, the highest neutron exposure of 0.44 mbarn^^ 
reached among the three representative layers of the ^^C pocket assumed in Gallino et al. 
(1998, see their Fig. 6), has to be considered as the maximum allowed for stars of solar 
metallicity when the mass fraction of ^^C in the He intershell Xq"2^''*^^" ^ 0.23. 

The recycling/leaking effect is affected by the choice of the ^^C(p,7)^'^N rate. When the 
lower rate from Caughlan & Fowler (1988) is used ^^C is less destroyed by proton captures. 
For example, at X^'^^ = 0.003, with the lower rate we obtain an increase of the value of r 
by 18%. 

As mentioned in the introduction, for AGB stars of metallicity lower than solar the 
s-proccss efficiency is higher for the same initial ^^C amount. This is because there arc fewer 
iron nuclei seeds per neutron and fewer light isotopes to act as neutron poisons. However, as 
well as ^^C, also the ^^N nuclei in the pocket are of primary origin; i.e. they do not depend on 
the stellar metallicity but are produced starting only from the H and He originally present in 
the star. Hence, the effect of the ^^N neutron poison will also be important in stars of lower 
metallicity and detailed computations are needed to address this point. Also rotation could 
affect the neutron exposure because in this case mixing continues after H-burning reignition 
and ^^N is mixed down into the ^^C-rich region during the interpulse period. This has the 
effect of strongly lowering the neutron exposure (Herwig, Langer, & Lugaro 2002). 

We conducted other simulations in the same 19*^^ interpulse period by artificially chang- 
ing the intershell amount of ^^C. As shown in Fig. 5 a nearly linear relationship exists 
between the final maximum value of r in the pocket and the intershell ^^C because the latter 
determines the ^^C and ^^N profiles and how many protons are recycled. The mass fraction 
of ^^C in the intershell varies from pulse to pulse (see e.g. Fig. 6 of Schonberner 1979, and 
Fig. 9 of Boothroyd & Sackmann 1988c) but it typically reaches a constant value around 
0.20 - 0.25 for the advanced pulses where the TDU occurs in the FRANEC and MSSSP 



-16- 



computations, while Xqw is only ~ 0.005. In the EVOL case instead, when hydrodynamic 
overshoot is included and intershell dredge up is at work, the intershell mass fraction of 
^^C varies from 0.30 to about 0.50 and then decreases towards an asymptotic value around 
0.40. The intershell mass fraction of ^^O increases to about 0.20 and then decreases towards 
about 0.15 (see Fig. 11 of Herwig 2000). The maximum neutron exposure in the pocket will 
then vary from pulse to pulse following the ^^C abundance. 

We have determined here that, given the reaction rates of proton and neutron captures 
on ^^C, ^"^C and ^^N , for a given metallicity and for a given intershell abundance of ^^C there 
is a maximum value of the neutron exposure achievable in the pocket, and hence a possible 
variation of neutron exposure is allowed to be in the range where r < Tmax- At solar 
metallicity the maximum s-process efficiency is comparable to the maximum value used by 
Busso et al. (2001) with FRANEC to explain the spread of neutron exposures required by 
spectroscopic observation. 



3.2. The choice of the ^^C pocket for the s- process computations 

To follow the ^^C pocket with the s-process code we divided the ^"^C-rich layer at 
the envelope-core interface into a number of zones of suitable mass and neutron expo- 
sure. The nucleosynthesis in each of these zones is computed individually until the end 
of the interpulse period. The resulting final abundances in the pocket are obtained as 
final ^ {^Xf'"'^M^one)/^M,one, wherc is the final abundance of the isotope i in 

each zone, M^one the mass of each zone, and the sum is done over all the zones. 

To specify the actual ^^C-pocket profile to run the s-process code with the structure 
computed by FRANEC and MSSSP, an initial H profile needs to be chosen which relates the 
Xp X-axis in Fig. 4 to the Lagrangian mass coordinate. The extent in mass of the pocket 
and of each of its zones has to be considered as a free parameter in these computations, thus 
leaving open the question of the specific shape of the H-profile and the mixing processes that 
lead to the partial mixing. We follow the approach of Gallino et al. (1998) and take the same 
profile that is plotted in their Fig. 1. The total mass for the pocket is 5 x 10^^ Mq, only the 
region of the pocket where the ^^N abundance is lower than the ^^C abundance, and hence 
the s-process is efficient, is considered and more mass is given to the zone characterized by 
the lower neutron exposures. The ^^C profile is kept the same for each cycle, however, as 
discussed by Gallino et al. (1998, see their Fig. 6), the neutron exposure slightly decreases 
with the interpulse number. This is because the initial composition of the s-processed ma- 
terial in the ^^C pocket is different due to the interplay of the different phases through the 
whole evolution. 
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For the EVOL case, we used the ^^C pocket self-consistently computed at the 5 inter- 
pulse period with the overshoot parameter / = 0.128. The initial amount of ^^C, ^^C and 
^^N, together with the resulting final neutron exposure as functions of the position in mass 
are plotted in Fig. 6. The total mass of the pocket is 3.7 x 10~^ Mq. As noted in the previous 
section, the higher the initial ^^C amount, the higher is the maximum r allowed (see Fig. 
5). Because of intershell dredge up the intershell mass fraction of ^^C found with the EVOL 
code, at this interpulse, is 0.48. The maximum of the neutron exposure is 0.84 mbarn"^, 
double the value obtained with the FRANEC and MSSSP codes, for which ~ 0.23. 

We note that this maximum neutron exposure is about 20% higher than the value of 0.70 
mbarn"^ obtained with the tests presented in the previous section and shown in Fig. 5 for 
the same ^^C abundance in the He intershell. This is due to the fact that the profiles of Fig. 
6 have been obtained using the -'^^C(p,7)"'^^N reaction rate from Caughlan & Fowler (1988), 
which, as already shown in the previous section, produces a small increase in the neutron 
exposure. 



3.3. The neutron flux in the ^^C pocket 

To compare the neutron fluxes in the interpulse period as obtained with the three 
different evolutionary sequences we plot in Fig. 7 the temperature from the three stellar 
evolution calculations for a typical interpulse period and the resulting neutron density as 
computed by the s-process code in the layer of the pocket characterized by the higher neutron 
density. For FRANEC and MSSSP the temperature shown corresponds to a point in mass 
at about 3x10""^ Mq below the H/He discontinuity left by the third dredge up (TDU), for 
the EVOL code it corresponds to the point of the maximum of ^^C in the pocket, which is 
located at about 3xl0~^ Mq below the H/He discontinuity left by the TDU. In all three cases 
the temperature reaches 10^ K before the onset of the next thermal pulse so that all ^^C is 
consumed. In the FRANEC, MSSSP and EVOL cases a typical interpulse period lasts for ~ 
40,000, ~50,000 and ~70,000 yr respectively (sec Fig. 3). The time after the beginning of 
the interpulse period at which the temperature reaches 0.8 xlO^ K is about 10,000 yr for the 
MSSSP and the FRANEC codes and 20,000 yr for EVOL. This is because the temperature 
in the MSSSP and FRANEC cases has been taken at a deeper, hence hotter, layer than in 
the EVOL case, however, it will not produce any important difference in the process of the 
neutron release. Since all ^^C is consumed, the total integrated neutron ffux depends on the 
initial amount of ^^C and not on the way it burns. The neutron flux lasts for about 20,000 
yr in the three cases, only being shifted in time because of the different temperature trends. 
The neutron density peak reaches 1.5, 0.9 and 2.3 xlO^ n cm~^ using the FRANEC, the 
MSSSP and the EVOL codes respectively. These variations are unimportant because the 
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neutron density is too low in any case to affect any branching on the s-process path. Note 
that in the temperature range of interest, T ~ 0.8 — 1 x 10^ K, the ^^C(Q;,n)^^0 rate we 
use is about 50—70% lower than the one adopted in the NACRE compilation (Angulo et al. 
1999). This difference would not make much change in the overall result. If ^^C burns faster 
the maximum neutron density will be somewhat higher, yet still unimportant in affecting 
the s abundances of branching- dependent isotopes. 

In conclusion, the important difference among the three computations regarding the 
neutron flux in the ^^C pocket is the maximum neutron exposure reached within the region. 
As illustrated in detail in the previous sections the maximum neutron exposure is doubled 
in the EVOL case with respect to the FRANEC and MSSSP cases because intershell dredge 
up in EVOL produces a double intershell abundance of ^^C. In the EVOL computation the 
s-element distribution will be pushed toward heavy s-process elements. 

3.4. The neutron flux in the convective regions by the ^^Ne neutron source 

We already mentioned in the introduction that ^^Ne is produced in the He intershell dur- 
ing the pulse driven convection zone (PDCZ) by the ^"^N(q;,7)^^F(/5"'"i/)^*0 (Q;,7)^^Ne chain 
that converts to ^^Ne all the ^^N significantly present in the ashes of the H burning and 
ingested in the PDCZ. The ^^N in the H ashes derives from the conversion of initial CNO 
nuclei hence it depends on the metallicity. However, because primary ^^C (and also ^^O, if 
intershell dredge up is at work) is dredged up from the He intershell in the envelope, the 
abundance of ^^N in the H ashes also contains a primary contribution that grows with the 
number and the strength of TDU episodes (Gallino et al. 1998). In Fig. 8 we show the 
contribution of this primary component to the ^^N ingested in the PDCZs as computed by 
the different codes. In our computations the primary ^'^N that is present in the ^^C pocket, 
as discussed in §3.1, does not contribute significantly to the total ^^N ingested by the PDCZ 
because the mass of the pocket is very small compared to the mass of the material processed 
by the H-burning shell (the typical ratio M pocket / is around 1/15). 

The s-process code activates the ^^Ne(Q;,n)^^Mg reaction when the bottom temperature 
of the PDCZ grows above 2.5x10^ K. We took the rate of the ^^Ne(Q;,n)^^Mg reaction from 
Kappeler et al. (1994), excluding the elusive resonance at 633 keV and using the lower limit 
of the strength of the 828 keV resonance. In the temperature range of interest, T ~ 2 — 4 x 10^ 
K, this rate is within a few percent of that adopted in the NACRE compilation (Angulo et 
al. 1999). while it is about 30% higher than a most recent determination by Jaeger et al. 
(2001). Using this latest estimate would have a small effect on the results presented here, 
only changing little the abundances of the isotopes connected to the branching points that 
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are very sensitive to the activation of the ^^Ne(Q;,n)^^Mg reaction in the PDCZ, such as ^^Zr. 

To compare the neutron fluxes during the thermal pulse the flrst important parameter 
to consider is the maximum temperature at the bottom of each pulse driven convection zone. 
As shown in Fig. 3 these temperatures when computed by the FRANEC and MSSSP codes 
are basically the same. A maximum temperature of ~ 3.05 x 10^ K is found in the last PDCZ. 
With the EVOL code higher temperatures are achieved as an effect of the hydrodynamic 
overshoot applied at the base of the PDCZ and the maximum reaches ~ 3.20 x 10^ K. 
Hence, the ^^Ne(Q;,n)^^Mg neutron source works more efficiently when the computation of 
the stellar evolution includes overshoot. Not only is the maximum temperature higher in this 
case, it is also high for a longer time and the overall distributions in time of temperature and 
neutron density in a PDCZ change. This is illustrated in Fig. 9, where we plot the bottom 
temperature trends for the last PDCZ computed by the three different stellar evolution 
calculations, together with the neutron densities computed by the s-process code on the 
basis of the three evolutionary codes. The temperature profile in the PDCZ has a narrow 
peak and hence so does the neutron density. For the FRANEC and MSSSP computations 
the neutron exposures by the ^^Ne neutron burst are very small, about a factor of 20 smaller 
than those in the ^^C pocket. Note that this can be different for AGB star models of higher 
mass or of metallicity lower than solar (Z < Zq/2), where the temperatures at the bottom 
of the PDCZs are higher (see e.g. Lugaro et al. 1999). However, the neutron densities reach 
up to 10^° n cm~^ and branching points are extremely sensitive to this neutron flux (Gallino 
et al. 1998). In the EVOL case the neutron exposure in the PDCZs reaches almost 0.1 
mbarn"^, which is a value high enough to influence the production of the s-process elements 
of the Sr peak. Also the neutron density is higher, up to about 6 xlO^^ n cm~^, hence we 
expect that in the EVOL computation neutron-rich isotopes produced by branchings on the 
s-process path will be produced in greater amounts. An important example is ^^Zr: the 
branching at the unstable ^^Zr is open when the neutron density is higher than ~ 5 x 10* n 
cm~^. From Fig. 9, we see that this occurs for about a year in the FRANEC and MSSSP 
cases and for about two years in the case of the EVOL code. The tail of the temperature 
distribution is also important in determining some flnal abundances of isotopes involved in 
branchings because of the freeze out effect (Cosner, Iben, & Truran 1980). In any case only 
a small amount of ^^Ne burns: for the PDCZs presented in Fig. 9, 2.6%, 1.1% and 8.9% of 
^^Ne burns in the FRANEC, MSSSP and EVOL cases, respectively. 
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3.5. The resulting s-process distributions 

In Fig. 10 wc show the production factors with respect to solar abundances from Anders 
& Grevesse (1989) of the s-only isotopes in the He intershell material averaged over the 
dredged up mass: x/™"' = Yl!j=i-^i-^^TDu/ Yl!j=i^TDU- ^^^'^ TDU episodes, 

Xl is the abundance of the isotope i in the He intershell after the corresponding pulse 
with TDU j and M^^jj is the mass dredged up. Also included are the production factors 
of ^^Sr, ^^^Ba and ^°^Pb, representative of neutron magic nuclei at N — 50, 82 and 126 
respectively, and of three neutron-rich nuclei most sensitive on branchings: ^^Kr, ^^Rb and 
^^Zr. The s-process abundances approach an asymptotic distribution after about 5 to 10 
TDU episodes. The production factors in the He intershell are up to a few hundreds and 
result in enhancements with respect to Fe of up to an order of magnitude in the envelope. 

It has been noticed in earlier studies, and pointed out again in recent work (see for 
example Busso et al. 1999), that the s-process abundance distribution is grossly determined 
once the neutron capture cross sections and the total number of neutrons available are 
given. In other words the neutron exposure determines the steps in the s-process abundance 
distribution at neutron magic nuclei and the slope of the distribution between the steps. 
This means that the general trends of the distributions shown in Fig. 10 depend on the 
choice of the neutron exposures in the ^^C pocket. Although the FRANEC and MSSSP final 
distributions are very similar since they use the same ^^C pocket, they show small differences: 
the more noticeable of these is that the production factors obtained with the MSSSP code 
are slightly lower than those obtained with the FRANEC code; moreover, the distributions 
are slightly different. This can be explained by the different structure inputs, mainly the 
mass of the convective regions, which is higher in the MSSSP case (see Fig. 3) resulting in 
a higher dilution of the material from the pocket in the He intershell. 

For the EVOL case the situation is different because of the higher neutron exposures in 
the ^^C pocket that are obtained owing to the large abundance of ^^C in the He intershell 
produced by intershell dredge up. While in the FRANEC and MSSSP cases the isotopes 
between the neutron magic nuclei ^^Sr and ^^^Ba are more enhanced than the isotopes 
between the neutron magic nuclei ^^^Ba and ^'^^Pb, the situation is reversed in the EVOL 
case because of the higher neutron exposures in the ^^C pocket. Also at the very end of 
the s-process path, ^°^Pb and ^°^Bi are produced in much greater amounts in the EVOL 
case. As noted in the previous section in the EVOL case the neutron exposure in the 
PDCZs reaches 0.1 mbarn"^, about five times more than the maximum value achieved with 
FRANEC and MSSSP. In this situation the abundances of isotopes with 70 ^ A < 100 are 
mostly produced by the activation of the ^^Ne neutron source. In fact, a simulation of the 
EVOL case performed by excluding the ^^Ne neutron source gave production factors of 2, 
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15, 30 for ™Ge, ^^Sr and ^^'^Ru, instead of 99, 96 and 50, respectively. Furthermore, the 

abundances of isotopes not affected by branching points from ^^^Ba to ^^'^Pb resulted to be 
lower by about 50%, and the production factor of 208p|^ ^^^g lowered by a factor of 3. The 
distribution in the region between 100 < A < 138, which is almost flat for FRANEC and 
MSSSP, is instead characterized in the EVOL case by a minimum at the s-only ^^^Sn. Such 
minimum is still present but much less accentuated in the simulation where the ^^Ne neutron 
source was excluded. 

The different modality of activation of the ^^Ne neutron source produces variations in 
the flnal results of isotopes involved in branchings. We include in Fig. 10 the production 
factors of three interesting isotopes produced by branchings on the s-process paths: ^^Kr, 
which is sensitive to the branching at ^^Kr; ^^Rb, which is sensitive to the branchings at 
^^Kr and ^^Rb; and ^^Zr, which is sensitive to the branching at ^^Zr. Because of the different 
neutron density profiles in the PDCZs (see Fig. 9), ^^Kr, ^^Rb and ^^Zr are produced in 
much higher amounts in the EVOL computation. The fact that the MSSSP case produces 
less *^Kr, ^^Rb and ^^Zr with respect to the FRANEC case is due to the fact that the MSSSP 
model has fewer PDCZs with the highest temperature than the FRANEC because of the 
effect of the higher mass loss rate adopted. 

Another interesting isotope is the s-only ^^^Gd whose production depends on the branch- 
ing at ^^^Sm. This isotope is actually positioned on the neutron-poor side of the branching 
and hence it may be surprising that is more produced in the EVOL case which has the 
higher neutron density. However, this is explained by the freeze out effect: the abundances 
of isotopes connected with branching points freeze out at the time when "aTrem < 1, where 
a is the largest cross section involved in feeding the s-only isotope characterizing the branch- 
ing, and Trem IS the remaining neutron exposure until the end of the pulse" (Kappeler et al. 
1990). Hence, for higher cr, such as in the case of the branching at ^^^Sm, final abundances 
are determined late during the neutron flux when the neutron density is low. The long tail 
of the neutron density distribution in the EVOL case (sec Fig. 9) allows enough time for 
the isotopes connected to the "high-o"" branchings such as ^^^Gd to determine most of their 
abundances at a late time during the neutron flux. 



4. Compcirison with observational constraints 

From the earliest (Ulrich 1973; Iben 1975b) to more recent studies (Arlandini et al. 
1999) a major constraint for the s process in AGB stars has been the reproduction of the 
solar abundance distribution of the s-only nuclei belonging to the main component (i.e. the 
s-process isotopic distribution for atomic mass number A > 90). From Fig. 10 it is clear 
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that our resulting abundance distributions are non-solar since the production factors of s- 
only isotopes are not all the same. It has been already pointed out by Gallino et al. (1998), 
making use of the FRANEC code, that AGB stars of solar metallicity with our choice of 
the ^^C pocket do not reproduce the main component. We can extend this conclusion to 
the MSSSP and EVOL cases. However, since the abundance distribution mostly depends on 
the neutron exposures achieved in the ^"^C pocket, which in turn depend on the ^^C amount 
(hence also on the intershell ^^C abundance) and, as previously discussed, on the stellar 
metaUicity, roughly following the rule r ~ ^^C/Z (Clayton 1988), we can obtain a best fit 
to the main component by adjusting these parameters. In the FRANEC and MSSSP cases 
we need to produce more of the heavier elements beyond the step at ^^^Ba by increasing the 
neutron exposure by about a factor of two. To accomplish this using the same ^^C pocket we 
have to decrease the stellar metallicity to half solar (Gallino et al. 1998). For the EVOL case, 
on the other hand, we would need to decrease the neutron exposures in the ^^C pocket. We 
could try to accomplish this by using a smaller overshoot parameter at the base of the PDCZ, 
hence having a lower ^^C amount in the intershell, or by including some rotationally induced 
smearing, which would mix ^^N in the ^^C pocket during the interpulse thus lowering the 
neutron exposure. In the EVOL case, though, we could still face great problems connected 
with the activation of the ^^Ne neutron source, such as the overproduction of ^^Kr, ^^Rb, 
96Zr, and ^^^Gd. 

However, we will not pursue this point because "it is a basic premise of the Burbidge 
et al. (1957) theory of nucleosynthesis in stars" (Clayton, Fowler, & Zimmerman 1961) that 
solar abundances are expected to come not from a single star, but rather to be built up 
by galactic chemical evolution. The more realistic approach to the origin of s elements in 
the solar system has been presented by Travaglio et al. (1999, 2001) and involves the study 
of the galactic chemical evolution of heavy elements making use of yields from AGB stars 
of different masses and metallicities. Making use of the stellar yields computed with the 
FRANEC and s-process code, these authors have shown that AGB stars of low mass can 
contribute satisfactorily to the Galactic abundances of the heavy s elements. The roles of the 
s-process weak component from massive stars (Lamb et al. 1977; Raiteri, Gallino, & Busso 
1992) and of neutron capture nucleosynthesis in supernovae of Type II have to be reviewed 
in order to fully understand how many realistic contributions to the light s elements up to 
A ~ 90 in the solar system are to be expected (Gallino et al. 2001). 

Hence we compare our model predictions with other observational constraints. Such 
comparison is not intended to be complete since we are only considering one stellar model 
of given mass, metallicity and ^^C pocket. We aim at understanding where the three dif- 
ferent evolutionary sequences stand with respect to each other when related to observa- 
tional constraints. The s-elements enhancement at the stellar surface for AGB stars of solar 
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metallicity are up to a factor of ten (Busso et al. 2001), which is typically matched by our 
models after dilution of the He intershell material with the convective envelope. The ratio 
[hs/ls]=log[(hs/ls)meas/(hs/ls)Q], where 'Is', which is the average of light s-process elements 
(Y and Zr), and 'hs', which is the average of heavy s-process elements (Ba, La, Nd and Sm) 
is a good indicator of the distribution of s-elements, and varies in the range —0.6 - for 
AGB stars of solar metallicity (Busso et al. 2001; Abia et al. 2002). The trends of our model 
predictions for [hs/ls] follow from what we observed in §3.5: in the FRANEC and MSSSP 
cases light s-process elements are more produced than heavy s-process elements and [hs/ls] 
is negative, while the opposite is true for EVOL. In this case, [hs/ls] is ~ 0. In principle 
stars with [hs/ls] ~ cannot be covered by the FRANEC and MSSSP predictions at solar 
metallicity, however, they can fit into these model at slightly lower metallicities, hence within 
observational uncertainties. 

Together with spectroscopic observations of AGB stars we also consider recent high- 
precision laboratory measurements of isotopic ratios of ,s-process elements in presolar silicon 
carbide grains (SiC, see review by Hoppe & Ott 1997). The vast majority (>90%) of SiC 
grains from the Murchison meteorite, the mainstream grains, most likely condensed in the 
mass losing envelopes of carbon stars. Note that the condition C/0>1 has to be satisfied for 
SiC molecules to condense. In particular, stars of low mass and around solar metallicity have 
been identified as the most likely parent stars of mainstream SiC (see Lugaro et al. 1999, for 
a thorough discussion). Measurements have been performed on aggregates of grains for Sr, 
Ba, Nd, Sm and noble gases. Recently, a new technique, developed at the Fermi Institute 
of Chicago (Nicolussi et al. 1997), has allowed measurements of isotopic anomalies of Sr, Zr, 
Mo and Ba in single SiC of the size of few /xm. The ratios ^^Sr/^^Sr and ^^^Ba/^^^Ba are 
sensitive to the neutron exposure. Measurements in aggregates of SiC grains indicates that, 
on average, the s-process distribution produced in SiC grain parent stars is accounted for by 
a non-solar distribution produced by neutron fiuxes in the pocket of about half the efficiency 
needed to reproduce the solar main component, such as that obtained in the FRANEC and 
MSSSP cases. The EVOL case produces instead too high ^^^Ba/^^^Ba and too low ^^Sr/^^Sr 
ratios. The ^^C pocket generated by the stellar model presented by Cristallo et al. (2001) 
has a neutron exposure profile with the same maximum s-process efficiency as we used in 
the FRANEC and MSSSP computations, as Xg^*!''^'^^" ~ 0.23, but with a distribution in 
mass more weighed towards the regions with higher s-process efficiency. This results in an 
s-process distribution slightly pushed towards heavy s elements, but little changes in the 
final results. The main difference is that ^^Sr/*^Sr ratio increases going outside of the range 
observed in SiC. 

All the above considerations, however, depend on the main model uncertainty: the 
^^C pocket. For example, as already mentioned, making use of the FRANEC code, Busso et 



-24- 



al. (2001) showed that one has to assume a spread in the ^^C-pocket s-process efficiency to 
cover the range of spectroscopic observations at any metalhcity. The same conclusion was 
reached by Lugaro et al. (2002) when considering single SiC grain data. It would be plausible 
to decrease the s-process efficiency in the pocket by means of stellar rotation as shown by 
detailed modeling (Herwig et al. 2002) so that the effect of rotation may naturally explain 
the needed spread of ^^C-pocket efficiencies at each metalhcity. This effect would also lower 
the neutron exposure in the pocket found by EVOL. 

Let us consider a few observational constraints that are sensitive to the high neutron 
exposure deriving from the the activation of the ^^Ne in the PDCZ. The high efficiency of the 
^^Ne neutron source produces in EVOL a situation similar to that occurring in intermediate- 
mass stars and in low-metallicity stars. The fact that overshoot at the base of the PDCZ 
induces a higher occurrence of a captures on ^^Ne is in contrast with spectroscopic observa- 
tions of Mg, which have shown that Mg isotopes are present in solar proportion in AGB stars 
of about solar metallicity (Clegg, Lambert, & Bell 1979; Smith & Lambert 1986). In the 
FRANEC and MSSSP cases we have a modest increase of the ^SMg/^^Mg and ^^Mg/^^Mg 
ratios in the envelope, from the solar values of 0.13 and 0.14 respectively, to about 0.16. In 
the EVOL case these ratios increase more abruptly, up to values around 0.45 (see also Herwig 
2000). The [Rb/Sr] ratio in AGB stars is typically negative (Lambert et al. 1995; Abia et 
al. 2001). In the FRANEC and MSSSP cases the [Rb/Sr] ratio is always negative, while in 
the EVOL case is always positive, at odd with observations. Large depletions of ^^Zr/^^Zr 
are seen in single SiC (Nicolussi et al. 1997) indicating that neutron densities in the stellar 
sources of these grains must have been low. These measurements are compatible with the 
FRANEC and MSSSP cases. The EVOL model produces instead only large positive values 
of the ^^Zr/^^Zr ratio. All these effects depend on the choice of the free parameter / at the 
base of the PDCZ. 



5. Conclusions 

We presented detailed computations of the s-process abundance distribution of heavy 
elements for a 3 Mq star of solar metallicity making use of three different evolutionary 
codes. One of them includes a time-dependent hydrodynamic overshoot mechanism and 
finds the ^^C pocket self-consistently. Hydrodynamic overshoot mechanism at the base of 
the convective envelope is a likely process to produce a ^^C pocket. The alternative mixing 
process related to rotation may limit the neutron flux by continued mixing of ^"^N into the 
s-process layer during the interpulse phase (Herwig et al. 2002). A combination of both 
processes should be investigated in the future. It is of course conceivable that other mixing 
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processes could be involved, which may have complementary or similar properties to those 
of the hydrodynamic overshoot mechanism considered here. Overshoot at the base of the 
convective envelope generates higher TDU efficiencies making it possible to create carbon 
stars at very low luminosities. The inclusion of overshoot at the base of the PDCZ, and the 
consequent intershell dredge up, help to explain the observed carbon and oxygen abundances 
of H-deficient central stars of planetary nebulae. 

The two computations that do not include overshoot do not self-consistently find any 
proton diffusion at the end of TDU and we assumed an identical ^'^C pocket for these cases. 
The final results of these two cases arc remarkably similar and we can conclude that the 
main uncertainty in s-process predictions is the ^^C pocket and not the intrinsic differences 
among different codes when no overshoot mechanism is included. In the case that includes 
hydrodynamic overshoot, because of penetration of the intershell convective region into the 
CO core, (i) the abundance of ^^C in the intershell and hence the maximum total neutron 
exposure achieved in the ^^C pocket are higher than in the other cases, and (ii) the neu- 
tron flux in the PDCZs is stronger. In summary, overshoot at the base of the convective 
envelope produces a ^'^C pocket if large efficiencies, i.e. values of the parameter /, are as- 
sumed; overshoot at the base of the PDCZ inffuences the activation of the ^^C(a,n)^^0 and 
^^Ne(Q;,n)^^Mg neutron sources producing a more efficient neutron flux in both situations. 
The final resulting abundance distributions of heavy elements reflect such features of the 
neutron fluxes: in the case with overshoot included, the heaviest elements, such as Sm, are 
produced in greater amount than lighter elements such as Sn. At the end of the s-process 
path ^°^Pb is produced in a higher amount, as are also isotopes fed by branching points, 
such as ^^Kr, ^^Rb and ^^Zr. 

With our choice of the ^^C-pocket proflle the results from the two cases without over- 
shoot give a good match to several observational constraints. The strong intershell dredge 
up connected with overshoot at the base of the PDCZs is in several situations in contrast 
with observational constraints. We still need to expand this study to a range of masses and 
metallicities. Within the EVOL case, we also need to extend our calculations in the space of 
the free overshoot parameter /, as applied at different boundaries and possibly also different 
times during the evolution, in order to test which range of / could provide us with the closest 
match to observational constraints. For example, a preliminary test with /=0.008 applied 
at the base of the PDCZ show a signiflcantly reduced efficiency in the activation of the ^^Ne 
neutron source. 
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FIGURE CAPTIONS 

Fig. 1. — Structure evolution through two thermal pulses of a 3 M© star model of solar 
metallicity computed with the MSSSP stellar evolutionary code (see §2). The y-axis is the 
Lagrangian mass coordinate in solar masses. Convective regions are white and it is possible 
to recognize the pulse driven convection zones during the the 19**^ and 20**^ pulse, as well 
as the bottom of the convective envelope. The grey lines in the radiative regions represent 
the mass shells. The region with the higher mesh density is the H-burning shell. The 
^^C pocket is a local phenomenon involving a mass region of the order of half a tick-mark 
(or less) on the y-axis. The x-axis is a representation of time by model number (from the 
Monash nucleosynthesis post-processing code, sec §3). The corresponding approximated time 
intervals are as follows: starting from t=0 at the end of the 19''^ PDCZ, the 10*'^ TDU lasts 
for about 160 yr, the formation of the ^^C pocket occurs within t=2500 yr, the ^^C(a,n)^^0 
reaction starts at t=ll,000 yr and by t=42,000 yr ^^C nuclei are depleted to a number lower 
than 10-^ The ingestion of the pocket in the PDCZ occurs at t=52,000 yr. The PDCZ lasts 
for about 40 yr, with about 4 yr during which the bottom temperature is above 2.5 xlO^ K 
and the ^^Ne(Q;,n)^^Mg reaction is marginally activated. 

Fig. 2. — TDU mass and envelope features as functions of the mass of the H-exhausted 
core for the 3 Mq and solar metaUicity models computed with the FRANEC (open squares) , 
MSSSP (open triangles) and EVOL (open circles) codes. Each point corresponds to a thermal 
pulse followed by TDU. These models were used to compute the s process. For each sequence 
the full symbol represents the starting point from which the evolution develops following the 
dotted lines connecting the open symbols. Note that the core mass always increases with 
time in the FRANEC and MSSSP cases, while in the EVOL case it decreases with time after 
6 pulses. For each model we show: a) the mass dredged up by each single TDU episode, 
b) the envelope mass, c) the cumulative mass dredged up, and d) the C/0 ratio in the 
envelope. 

Fig. 3. — Structural and thermodynamical characteristics as function of the core mass. 
Each point corresponds to a thermal pulse followed by TDU. For each sequence the full 
symbol represents the starting point from which the evolution develops following the dotted 
lines connecting the open symbols, a) The interpulse period, note that in the EVOL case 
the interpulse period increases with the pulse number, while for FRANEC and MSSSP it 
decreases, b) The mass of the PDCZ, which decreases with time in all cases, c) The overlap 
factor r between subsequent PDCZs, which decreases with time in all cases, and d) the 
variation of the maximum temperature at the bottom of the PDCZs, which increases with 
time in all cases. 
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Fig. 4. — Resulting mass fractions of ^^C, ^^C and ^^N as functions of the initial mass fraction 
of protons introduced below the H/He discontinuity during the simulations performed with 
the Monash nucleosynthesis code. The protons are introduced after the 10*'^ TDU episode 
of the 3 Mq star of solar metallicity computed with the MSSSP code. Also plotted is the 
corresponding neutron exposure at the end of the interpulse period after all the ^^C has 
burnt. 

Fig. 5. — Maximum neutron exposure in the pocket as a function of the He intershell ^^C from 
the simulations performed at the lO*'^ TDU episode of the 3 Mq model of solar metallicity 
computed with the MSSSP code. The full hexagon indicate the true Xq^^^'^^^^^ value at such 
interpulse resulting with MSSSP. 

Fig. 6. — Abundance profiles in the ^^C pocket found self-consistently by the EVOL code 
at the 5**^ interpulse period with the overshoot parameter / = 0.128. The large tick-marks 
on the X-axis measure 10~^ Mq intervals. Also shown is the resulting r at the end of the 
interpulse as subsequently computed by representing the abundance profiles with ten zones 
and performing one-zone computations with Monash nucleosynthesis code. 

Fig. 7. — a) Temperature from the three evolutionary codes, and b) neutron density as 
computed by the s-process code in the zone 3 of the ^^C pocket with the higher neutron 
density. The zero points in time represent the time from the start of the interpulse period, 
about 10,000 yr, when in the MSSSP case T = 0.79 X 10*^ K. Matter densities also grow with 
time from 2 x 10^ g cm~^ to 10*^ g cm~^ for FRANEC and EVOL, and to ~ 1.5 x 10"^ g 
cm-3 for MSSSR 

Fig. 8. — Primary contribution to the ^^N in the H-burning ashes as a function of thermal 
pulse number. Primary ^^N comes from the dredge up of ^^C (and ^^O in the EVOL case) in 
the envelope followed by H-shell burning. To compute the total mass fraction of ^'^N engulfed 
by each PDCZ these values have to be added to the constant value of 0.0128 adopted in the 
s-process code for the H-burning ashes of solar metallicity. 

Fig. 9. — a) The temperatures at the bottom of the PDCZ and b) the neutron densities as 
functions of time during the last thermal pulse computed in the three cases: the 25*^ with 
TDU for FRANEC (sohd lines), the 19*^ with TDU for MSSSP (long-dashed fines), and tfie 
13*'' witfi TDU for tfie EVOL code (sfiort-dasfied lines). Tfie zero in time corresponds to tfie 
instant in wfiicfi tfie bottom temperature in tfie convective sfiell reacfies 2.5x10^ K. Bottom 
densities decrease exponentially witfi time from ~10^ to ~10^ g cm~^. Tfie variation of 
temperature, matter density and neutron density witfi pulse number for tfie FRANEC case 
can be found in Figs. 10, 11 and 13 of Gallino et al. (1998). 
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Fig. 10. — Distribution of production factors of s-only isotopes (open symbols) in the He 
intershell material averaged over all TDU episodes and weighed by the mass dredged up 
(see text). The production factors of the neutron-rich isotopes *^Kr, ^''Rb and ^^Zr, which 
arc affected by branchings, are also shown as full symbols. The vertical lines are eye guides 
drawn at the stable neutron magic nuclei on the s-process path: ^^Sr, ^^^Ba and ^''^Pb, whose 
production factors are also plotted using open dash-lined symbols. The production factor of 
the s-only isotope ^^^Gd is indicated by an arrow (see text for discussion). 
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